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Abstract. The Garching-Bonn Deep Survey (GaBoDS) is a virtual 12 square degree cosmic shear and cluster 
lensing survey, conducted with the WFI@2.2m MPG/ESO telescope at La Silla. It consists of shallow, medium 
and deep random fields taken in _R-band in subarcsecond seeing conditions at high galactic latitude. A substantial 
amount of the data was taken from the ESO archive, by means of a dedicated ASTRO VIRTEL program. 
In the present work we describe the main characteristics and scientific goals of GaBoDS. Our strategy for mining 
( the ESO data archive is introduced, and we comment on the Wide Field Imager data reduction as well. In the 

second half of the paper we report on clusters of galaxies found in the background of NGC 300, a random archival 
field. We use weak gravitational lensing and the red cluster sequence method for the selection of these objects. 
Two of the clusters found were previously known and already confirmed by spectroscopy. Based on the available 
data we show that there is significant evidence for substructure in one of the clusters, and an increasing fraction of 
blue galaxies towards larger cluster radii. Two other mass peaks detected by our weak lensing technique coincide 
with red clumps of galaxies. We estimate their redshifts and masses. 
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1. Introduction The amplitude of the distortions induced into galaxy 

images is of the order of a few percent for lensing by large- 
Light bundles from distant background sources are dis- scale structure, and up to ten percent in the case of clus- 
torted by the tidal gravitational fields of the interven- ter lensing, depending on the mass of the cluster under 
ing matter distributed along the line of sight. This shear consideration. The intrinsic ellipticities of galaxy images 
can be statistically measured by determining the image dominate the noise of the shear signal. Thus, a weak lens- 
ellipticities of those background sources. Since it is lin- ing analysis requires deep data with good PSF properties 
early related to the underlying matter distribution, the in order to measure enough galaxy ellipticities with sum- 
shear can be used to measure statistical properties of the dent accuracy. A large field of view and independent lines 
large-scale structure in the Universe (cosmic shear). It f sight guarantee a representative view of the universe, 
can also be used to detect clusters of galaxies, without beating down the effect of cosmic variance. Deep wide field 
assuming any hydrostatic equilibrium, symmetry or rela- surveys provide this kind of data. 

tion between luminous and dark matter. For recent re - The Garching-Bonn Deep Survey (hereafter: GaBoDS) 

vie ws about weak gravitational lensi ng see|Melliei| (|l99g|), is a mostly virtual survey and covers a sky area of more 

and|Bartelmann fc Schneider! (|200l|). than 12 square degrees with the Wide Field Imager 1 at 

the MPG/ESO 2.2m telescope. Compared to other tele- 



Send offprint requests to: M. Schirmer 1 Hereafter we will use the term 'WFF for any Wide Field 

* Based on observations made with ESO Telescopes at the Imager instrument, and 'WFIQ2.2' when refering to the Wide 
La Silla Observatory Field Imager at the MPG/ESO 2.2m telescope. 
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scopes (VLT, CFHT, CTIO, WHT, KPNO) that have 



(Bacon et all 12000. 


2002; iMaohetall 2001; 


Kaiser ct al. 


2000; 


Ivan Waerbeke et al. 




2000, 


l200ll 




Wittman et al. 



20001) . the 2.2m telescope is small in size, yet it offers an 



outstanding image quality. It allows for an independent 
confirmation of previously obtained results, and the de- 
tection of possibly remaining systematics in the data as 
well as in the measurement process. 

In this paper we use weak gravitational lensing to de- 
tect dark matter haloes by looking for co herent shear pat- 
terns in galaxy images llSchneideiill99^) . In this manner 
galaxy clusters are detected directly relying on their most 
fundamental property, their mass, which is independent 
of their dynamical state and luminosity. Since structure 
formation in the Universe is highly sensitive to the cos- 
mological model, shear-selected samples of dark matter 
haloes will yi eld new insights into the process of cluster 
formation fseelErben et all 120001 iMargoniner et alll200.1 ; 



Mivazaki et all l2002t iDahle et all l2002t IWittman et al. 
200 ll l2002i for shear-selected hal oes). Based on simula- 
tions of the large-scale structure l)jain fc van Waerbekel 
2000) as well as on prediction s from Press-Schechter mod- 
els l)Kruse fe Schneider! |2000|) one expects some 10 shear- 
selected dark matter haloes per square degree. Therefore, 
the number of mass peaks in the GaBoDS data is expected 
to be of the order of 100. The mass of these peaks strongly 
depends on their redshift, as is shown in Sect. 5. 

This work is organized as follows. In Sect. 2 we discuss 
our strategy for data mining the ESO archive. An outline 
of the main characteristics of the GaBoDS fields follows in 
Sect. 3, together with a compact overview of our WFI data 
reduction process in Sect. 4. Subsequently we report on the 
detection of galaxy clusters by me ans of weak lensing and 
the red cluster sequence method ( Gladder fc YeeL [2OOO) 
in the NGC 300 GaBoDS field, and give a more detailed 
analysis for one of the clusters therein. We conclude in 
Sect. 6. 

2. Data mining the ESO archive 

2.1. Our ASTROVIRTEL program 

For the above mentioned weak lensing analysis we aimed 
for a sky coverage of at least 10 square degrees for 
GaBoDS. Only about 3 square degrees were finally ob- 
served, however, in 20 allocated nights of our own GO pro- 
gram due to unfavourable weather conditions. Including 
available data from the EIS 2 Deep Public Survey and 
COMBO-17 3 left us with about 4 missing square degrees. 
A manual search in the large ESO archive turned out to 
be unfeasible, since the search engine available at that 
time did not allow for filtering the data with respect to 
our requirements. The only usable fields we knew before- 
hand in the archive were the ones from the Deep Public 
Survey and one pointing of the Capodimonte Deep Field. 

2 ESO Imaging Survey 

3 MPIA Heidelberg 



Other WFI@2.2 data such as the five COMBO-17 fields 
were taken during MPG time before the 2.2m telescope 
started into service mode operation, and were thus not 
publicly available through the archive. Besides, the very 
low number of publications based on data taken with 
WFIQ2.2 did not allow for a direct identification of further 
usable data. We therefore proposed an ASTROVIRTEL 
program 5 6 , aiming at an en hancement of the querator 7 



search engine's capabilities I Pierfedericil l20o"il) . In order 
for a field to be included in the GaBoDS, the following 
requirements addressed by querator had to be met: 

— minimum exposure time in i?-band: ~ 5 kscc, 

— image seeing < l'.'O, 

— random, at high galactic latitude, 

— empty, i.e. avoidance of known very massive structures 
in or next to the field, no bright stars or large fore- 
ground objects inside the field. 

The first item in this list guarantees a high enough number 
density (> 10 arcmin" 2 ) of galaxies with securely mea- 
surable shapes. Furthermore, exposures must be taken in 
excellent seeing conditions, since the S/N for shape mea- 
surements decreases with the second power of the size of 
the PSF. The random character of the fields is needed to 
avoid a biasing towards certain types of objects, such as 
quasars or clusters of galaxies. In addition, it is ensured 
in this way that the fields sparsely sample the Universe 
along independent lines of sight, thus keeping the effect of 
cosmic variance small. The last point in the above item 
list is to further guarantee that a given pointing is usable 
for our purposes: massive structures such as large clus- 
ters of galaxies bias the search for unknown dark matter 
haloes as well as the measurement of a cosmic shear signal, 
which is about one order of magnitude smaller than the 
weak lensing signal of a large cluster. We are aware of the 
fact that this introduces a bias towards lower density lines 
of sight. By avoiding bright and large foreground objects 
we keep the usable area of an image as large as possible. 

2.2. About data quality 

It is clear that 'data quality' is a very ambiguous term, 
highly dependent on the science which is to be drawn from 
the data. Given the involved and time consuming reduc- 
tion of WFI data, we wanted to assess as best as possible 
the quality of archival data before any request or data re- 
duction. This is straight forward for items such as the total 
exposure time, available filters, presence of bright objects, 
availability of calibration frames, and the ambient condi- 
tions during which the observation were performed (moon; 
clouds; seeing, according to the seeing monitor). However, 



4 Only about a dozen papers based on WFI@2.2 data have 
been published up to April 2002. 

5 ASTROVIRTEL cycle 2: Erben et al., Gravitational lens- 
ing studies in randomly distributed, high galactic latitude fields 



http:/ /www.stecf.org/astrovirtel 
7 http:/ /archive. eso.org/querator 
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Table 1. Fields identified by our ASTRO VIRTEL pro- 
gram. The second column gives the exposure time that 
went into the coaddition (images with very bad PSF or 
seeing were rejected). Exposure times in parentheses indi- 
cate the total exposure time per association as identified 
by querator. 



Target 


Exposure time 


Image seeing 


B08p3 


4900 


^5600) 


0'.'88 


B08p2 


4900 


;5600) 


0'.'88 


B08pf 


4500 


;5600) 


0'.'88 


B0800 


7200 


;8300) 


0'.'88 


B08mf 


4500 


'6000) 


0'.'88 


B08m2 


4000 


;5800) 


0'.'88 


B08m3 


5400 


;6ooo) 


0'.'96 


Pal3 


5000 


'6120) 


l'.'O 


AMI 


7500 


;7620) 


l'.'O 


Comparison f 


5300 


;9300) 


0?97 


C04pf 


4000 


^5600) 


0'.'88 


C04p3 


4000 


;5600) 


0'.'83 


C04m2 


4000 


;4800) 


0'.'85 


C04m3 


4000 


;4800) 


0'.'85 


NGC300 


15100 


(« 25000) 


1'.'06 


C04p2 


4000 


;5600) 


0'.'86 


C04mf 


4000 


;4800) 


Or 86 


C04m4 


4000 


;4800) 


0'.'86 


C0400 


4800 


;5600) 





not all of those can be expressed in terms of numbers; 
some judgement has to be done upon visual inspection of 
the data. Other data quality issues, such as image see- 
ing or PSF properties can not be evaluated without more 
complex operations on the data itself. To address most 
of these points, the following demands were defined to 
ASTRO VIRTEL and querator: 

— measure the image seeing for all archived WFI@2.2 
data, 

— provide a preview facility, allowing quick visual inspec- 
tion of the data ahead of a request, 

— make available the according proposal abstract ('why 
was this particular observation done?') 

All but the last item were implemented. The image see- 
ing, crucial for a weak lensing analysis, was determined by 
automatically extracting all non-saturated stars from one 
of the eight WFI@2.2 chips, and averaging their FWHM. 
It was superior to the DIMM 8 seeing, since the DIMM on 
La Silla could not pick up local effects such as dome see- 
ing, zenith distance, flexure and temperature of the tele- 
scope's Serrurier structure, focus and tracking. On average 
we found the DIMM seeing to be O'.'l - 0'.'2 better than 
the image seeing for WFI@2.2. It was only in rare cases 
that the difference between the two drops below O'.'l. 

In this way a list of useful candidate fields was ex- 
tracted from the archive, minimizing the amount of un- 
usable data slipping into the reduction process. Fields 

8 Difference Image Motion Monitor 
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Fig. 1. Sky distribution of the GaBoDS fields. The size of 
the symbols depicts the covered sky area. All fields are at 
high galactic latitude. 

that were rejected at this late step suffered from scat- 
tered light or had very bad PSF properties 9 . Checking for 
PSF anisotropics in an image was very time consuming 
and thus not blindly performed on all data in the archive. 

Using the enhanced querator we found about 5 square 
degrees (20 pointings) of data in the ESO archive which 
passed our criteria, not counting the already known fields 
such as the ones from the EIS Deep Public Survey. This 
was twice the area we expected, but it must be noted 
that 15 of the pointings were done by a single observer, 
searching for trans-neptunian objects. The efficiency of 
our archive search with querator was about 75%, i.e. three 
out of four candidate fields were usable. Table lists the 
fields found together with their image seeing. 

3. Main characteristics of the GaBoDS fields 

Since we made heavy use of archival data, we were able 
to collect more than 12 square degrees of high quality R- 
band images with only 10 clear nights of own observations. 
For nearly half of the fields multicolour information is also 
available. As can be seen from Figs.^ an d[21 the GaBoDS 
fields can be split into a shallow part (6 square degrees, 
4-7 ksec total exposure time), a medium deep part (4.25 
square degrees, 8-11 ksec) and a deep part (2.0 square 
degrees, 13-56 ksec). The survey fields are randomly dis- 
tributed at high galactic latitude in the southern sky. Data 
in GaBoDS that was not taken by ourselves is summarised 
in Table H 

The image seeing in the stacked images is equal to or 
better than l'.'O, and our astrometric solution is accurate 
enough not to introduce artificial PSF anisotropies (see 
Sect. 4 for more details). The typical PSF quality of a 
stacked set of WFI@2.2 exposures in our survey can be 
seen in the lower right panel of Fig. [3J There we show 
the PSF anisotropies of a coadded image consisting of 57 
exposures, giving a total integration time of 27.1 ksec. The 

9 The anisotropy threshold depends on the overall quality of 
a particular data set. In general we rejected exposures with 
anisotropies larger than ~ 6%. 
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Table 2. GaBoDS data that was not taken by ourselves. 



2x10" 4x10 4 
Exposure Time [s] 



6x10" 



Fig. 2. Exposure times in GaBoDS. The peak at 56 ksec 
represents the Chandra Deep Field South (CDF-S). 

image seeing of this particular field is C/.'8, thus possible 
shortcomings in the astromctric algorithm could easily be 
seen. The rms PSF anisotropy merely amounts to 0.8%, 
making WFIQ2.2 very well suitable for weak gravitational 
lensing studies. The larger anisotropies in the field corners, 
especially in the lower left, are due to slightly tilted CCDs 
with respect to the focal plane (K. Meisenhcimer, private 
communication). Optical abberations play only a minor 
role for this instrument. 

Such quality of the data can only be achieved with a 
very carefully and frequently refocused telescope. The ef- 
fect of a slightly defocused telescope on the PSF is shown 
in the remaining three panels of Fig. [3J One can see that 
anisotropies become significant once the detector is out of 
focus. Furthermore, the PSF rotates by 90 degrees when 
one passes through the focal plane. This is characteris- 
tic for tangential and sagittal astigmatism. Still the PSF 
of WFI@2.2 is very homogeneous over the field of view, 
even when crossing chip borders. Thus larger dither pat- 
terns can be used for the observations, and a single smooth 
model can be fitted to the PSF in the stacked mosaic (see 
also Fig.|H|in Sect. 5). 

4. The art of WFI data reduction 

4.1. The GaBoDS pipeline 

The advent of multichip CCD cameras imposes new, high 
demands on data reduction. Pre-processing steps such as 
debiasing or flatfielding can be done independently on 
a chip-by-chip basis, allowing for efficient parallel pro- 
cessing on a multi-processor machine with sufficient disk 
space. Whereas these steps can be tackled using the same 



Sky coverage 


Data source, depth of field 


1.75 deg 2 


ESO Imaging Survey + GOODS 




(medium, one deep) 


5.00 deg 2 


ASTRO VIRTEL (mostly shallow, one deep) 


1.25 deg 2 


COMBO-17 survey (deep) 


1.25 deg 2 


ESO distant cluster survey 




(shallow, medium) 


0.25 deg 2 


NTT Deep Field, from the IR group at the 




MPI fur extraterrestrische Physik (deep) 


0.25 deg 2 


Capodimonte Deep Field (deep) 



well-known algorithms as for single chip cameras, an ac- 
curate astrometric and photometric calibration of WFI 
data requires techniques going well beyond those rou- 
tines. Different sensitivities of the CCDs and gaps between 
them lead to a very inhomogeneous exposure time and ac- 
cordingly noise in the coaddition. An accurate weighting 
scheme is essential in order to retain control over these 
effects in the stacked image. In the following we describe 
our approach to WFI data reduction. 

An almost fully automatic pipeline for WFI reduction 
was developed, based on existing software modules wher- 
ever possible, such as EIS drizzle 10 , the Terapix 11 soft- 
ware suite, FLIPS 12 (J.-C. Cuillandre, not yet publicly 
available), Eclipse 13 , Imcat 14 (N. Kaiser) and LDAC 15 
(E. Deul, Leiden Data Analysis Center) (catalogue for- 
mat and handling). IRAF was not used in the pipeline, 
apart from the drizzle coaddition, since it did not allow 
for efficient scripting and reduction of this specific kind 
of data. A number of bash shell scripts were wrapped 
around those tools (mostly stand-alone C programs), al- 
lowing for an efficient, flexible and almost fully automa- 
tized end-to-end reduction of WFI data in parallel mode. 
The usage of the pipeline is not restricted to the WFI@2.2, 
but data from other instruments such as FORSl/2@VLT, 
ISAAC@VLT, SUPRIMECAM@SUBARU, MOSAIC- 
I@CTIO, MOSAIC-II@KPNO and WFI@AAO has al- 
ready been successfully reduced. Supported architectures 
are Solaris, AIX, Linux and Dec-Alpha. The package will 
be released together with a detailed technical description 
(Erben et al., in preparation). 

Our pipeline was designed with the GaBoDS data in 
mind, i.e. empty fields at high galactic latitude, with a 
fairly large dither pattern of up to 3'.0, and a very large 
number of single exposures per pointing. Still, the usage of 
the GaBoDS pipeline is not resticted to empty fields only. 
This is shown in Sect. 5, where we present an analysis of 
galaxy clusters behind NGC 300. 



http:/ /www. eso.org/science/eis 
http:/ /terapix. iap.fr 

http:/ /www. cfht.hawaii.edu/~jcc/Flips/flips. html 
http: / / www.eso.org/projects/aot /eclipse 
http: / / www.ifa.hawaii.edu / ~ kaiser / imcat 
ftp: / / ftp.strw.leidenuniv.nl / pub /ldac /software 
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Fig. 3. PSF anisotropies for an intrafocal (upper left), focal (upper right) and extrafocal (lower left) exposure. The 
chosen scale for the stick length is the same for those three plots in order to show the increase in the anisotropies with 
respect to the focused exposure. The mean stellar ellipticities are 6.6%, 0.9% and 5.9%, respectively. The lower right 
panel depicts typical PSF anisotropies of a stacked WFI@2.2 i?-band image (~ 50 exposures with ~ 500 sec exposure 
time each). Note that the largest PSF anisotropy in the stacked image is as small as 0.8%. Compared to the other 
three PSF plots a different scale for the stick length was used in order to clearly show the anisotropies. 



4.2. The pre-processing 



Overscan correction, debiasing, fiat fielding: Apart 
from the astrometric, photometric and coaddition pro- 
cesses, all chips are processed individually, allowing for 
an easy parallelization of the code. 



During pre-reduction, any instrumental signatures 
present in the data are removed. This includes over- 
scan correction, bias subtraction and flat fielding with 
skyflats 16 . For the master biases all bias frames are median 



16 For WFl@2.2m telescope dome flats are inferior compared 
to skyflats. 
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combined with outlier rejection. Flat fields are combined 
in the same way, but each flat exposure is normalised to 
1 before the combination. Thus the different gains in the 
science images are still present after the flat fielding step. 
This is because we found that the chip-to-chip gain varia- 
tions can be better determined from a superflat, which is 
drawn from already flatfielded exposures. 

Creation of a superflat: Residuals of around 3-4% 
after normal flat fielding are common for WFI@2.2m, de- 
pending on sky brightness and the filter in use. A superflat 
is computed for correction of this effect, by median com- 
bining all images from a given observing run, using outlier 
rejection. Pixels that are a ffected by stars or g alaxie s are 
detected with SExtractor (B ertin fc Arnoutsl 119961) and 
masked beforehand. In doing so we prevent the bright ex- 
tended haloes around stars contributing to the superflat. 
The image constructed this way is then heavily smoothed, 
yielding an illumination correction image for every chip. 
All images are divided by their individual illumination cor- 
rection. They are also normalised to the same (the high- 
est) gain, which is accurately determined by comparing 
the modes of the individual superflats. Remaining residu- 
als in the sky background are typically below 2%, and in 
the case of absence of bright stars, even below 1%. 

Defringing: Besides the illumination correction a 
fringing model is calculated by subtracting the illumina- 
tion correction from the previous superflat. Hereby it is as- 
sumed 17 that the amplitudes of the observed fringes scale 
directly with the sky background. The fringing model is 
then individually scaled for each image and subtracted. 
Fringing in the i?-band is of the order of a few percent for 
WFI@2.2. Its correction makes sense only if more than 
~ 10 images were used in the calculation of the model. 
Otherwise the pixel-to-pixel noise in the fringing model is 
larger than the fringing amplitude itself. This would intro- 
duce more noise into the individual defringed images than 
what is taken out by the correction of the lower frequency 
fringes. Since we are interested in measuring shapes of 
faint and small galaxies we want to avoid additional pixel- 
to-pixel noise. However, all images in the GaBoDS survey 
were constructed from many more than just 10 images, 
thus the contribution of additional high frequency noise 
is small, and we profit from taking out the gentle fringing 
pattern. In the case of the redder /-band the fringing can 
be much more prominent, and is in general more difficult 
to remove. 

4.3. Astrometric calibration 

After the pre-processing a global astrometric solution and 
a global relative photometric solution is calculated for all 
images. This is where the reduction of WFI data becomes 



17 This assumption only holds for photometric nights. In non- 
photometric nights it can be impossible to remove fringing in 
the /-band, since the model is changing within minutes on 
scales smaller than the field of view of a single CCD. 



much more complicated than the one for single chip cam- 
eras. 

In the first step, high S/N objects in each image are 
detected by SExtractor: and a catalogue of non-saturated 
stars is generated. Based on a comparison with the USNO- 
A2 astrometric reference catalogue, a zero-order, single 
shift astrometric solution is calculated for each chip in ev- 
ery exposure. For a single-chip camera with a small field of 
view such an approach is often sufficient, but it no longer 
holds for multiple chip cameras with a large field of view. 
CCDs can be rotated with respect to each other, tilted 
against the focal plane, and in general cover areas at a dis- 
tance from the optical axis, where field distortions become 
prominent 18 . Fig.0]shows the difference between a zero or- 
der (single shift with respect to a reference catalogue) and 
a full astrometric second order solution per CCD. From 
this figure it is obvious that the simple shift-and-add ap- 
proach will not work for the whole mosaic. The issue is 
further complicated by the gaps between the CCDs and 
large dither patterns that are used to cover them. Thus, 
chips with very different distortions overlap. In addition, 
due to the large field of view, the observed patch of the 
sky must no longer be treated as a flat plane, but as a 
spherical curved surface. 

In the second step we use Mario Radovich's 
Astrometrix 19 (Terapix) package to fit third order polyno- 
mials to every chip in every exposure, in order to correct 
for the above mentioned effects and to find a global as- 
trometric solution. For this purpose all high S/N objects 
(stars and galaxies) detected in the first step are iden- 
tified with each other, including those from the overlap 
regions. The latter ones are most important in establish- 
ing a global astrometric (and photometric) solution, since 
the accuracy of available reference catalogues such as the 
USNO-A2 (0'.'2 rms) is insufficient for sub-pixel registra- 
tion. Thus the astrometric solution is determined from the 
data itself. The USNO-A2 is used only to fix the solution 
with respect to absolute sky coordinates within 0'.'2 rms. 
With Astrometrix we consistently achieve an internal as- 
trometric accuracy of 1 /20 - l/10th of a pixel (0'.'02-0'.'01) , 
thus the final PSF is mostly determined by the intrinsic 
PSFs of the single exposures (see Fig.[3J). Additional, arti- 
ficial seeing and PSF anisotropies are introduced into the 
stacked image on a very low level only, even for very large 
data sets such as the CDF-S, consisting of 150 WFI@2.2 
R-band exposures. This is a crucial requirement for our 
weak lensing analysis. 



A ZEMAX analysis (Philipp Keller, private communica- 
tion) for the 2.2m Ritchey Chretien telescope design and 
WFI@2.2's focal reducer shows that the radial field distortion 
for this layout increases with 8 — air 2 + a2r 6 . However, the 
total amplitude of this distortion is very small (~ 30 pixels). 
19 http:/ /www. na.astro.it/~radovich/WIFIX/astrom.ps 



M. Schirmer et al.: GaBoDS: Anatomy of galaxy clusters in the background of NGC 300 



7 



4.4. Photometric calibration 

Once an astrometric solution is found, a relative photo- 
metric solution is straight forward. Relative fluxes of ob- 
jects in different exposures and overlap regions are com- 
pared, allowing the calculation of relative photometric ze- 
ropoints for every chip and every exposure. Given two 
overlapping chips k and j, consider all i = 1...N objects 
and calculate the mean deviation of magnitudes K and J 

ZjWjjKj - Jj) 

Mk ' j = ' (1) 

with Wi = {<j\ + cTj) -1 , where a are the measurement 
errors of the corresponding magnitudes. Objects deviat- 
ing in Ki — Ji more than a user defined threshold are 
rejected. The relative zeropoints ZP\ for all N avaT over- 
lapping CCDs are determined by x 2 minimization with 
respect to ZP k : 

JVovcr 

X 2 = £ [M kJ - (ZP k - ZPj)} 2 . (2) 

k,j 

Finally, the relative zeropoints of photometric images are 
normalised so that their mean is zero. This approach as- 
sumes that the relative zeropoints are constant for every 
CCD 20 . An automatic absolute photometric solution is 
not yet implemented in the pipeline. 

4.5. A statistically optimized weighting scheme 

The effective exposure time for a stacked WFI mosaic is 
highly non-uniform. Read noise and the flatfields are chip- 
dependent, and gaps between the CCDs contribute further 
to the inhomogeneous depth of a stacked image. Applying 
a statistically optimized weighting scheme to the expo- 
sures during coaddition allows for a significantly improved 
object S/N ratio (up to a factor of ~ 1.5). 

We now describe our approach to assign an individual 
weight map for every science image. In a first step, a pixel 
is assigned its normalised skyflat value as a weight, which 
contains the information about relative gains between the 
CCDs and pixel-to-pixel sensitivities. Contrary to other 
methods, we do not detect 'bad' pixels (hot or cold pixels, 
pixels affected by cosmics, reflections or satellite tracks) 
by intercomparing all images in the stack, but on the indi- 
vidual images themselves. For the detection of permanent 
image defects, such as hot or dead pixels and bad columns, 
we use dark frames and superflat images. Affected pixels 
are set to zero in the corresponding weight map. Thus ev- 
ery chip in the WFI mosaic has its own basic weight map 
after the first step. Weight maps created in this way are 
the same for all exposures in a data set unless they were 
flatficlded with different skyflats or taken several weeks 
apart. 

20 Zeropoint variations for images taken with WFI@2.2 that 
were not superflat corrected are described in: 

:/ /www. Is. eso.org/lasilla/sciops/wfi/zeropoints . 
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Fig. 4. Difference in object position between a single-shift 
approach and a full two-dimensional second order astro- 
metric solution for the WFI@2.2. In other words, shown 
are the higher order terms needed for matching the CCDs 
to the sky. The patterns belonging to the left two chips are 
due to a rotation with respect to the mosaic. The maxi- 
mum position difference in the plot is about six pixels, still 
a fairly small value compared to other telescope designs. It 
becomes clear that a single, global distortion polynomial 
for all CCDs does not work. Instead, every CCD has to 
be treated individually. 

In a second step these weight images are adjusted in- 
dividually for every image. Remaining hot pixels and cos- 
mics are easily identified with SExtractor in conjunction 
with Eye (Terapix), since they appear much sharper than 
the stellar PSF even under good seeing conditions. Bright 
reflections and satellite tracks, however, need to be masked 
by hand, the only step in the pipeline which is not yet 
automatized. Moving objects like asteroids go unmasked 
and show up as dashed lines in the stacked image. During 
coaddition the individual weight maps are scaled with cor- 
rection factors for airmass and varying photometric con- 
ditions. Changing seeing conditions from image to image 
can be included on an optional basis, too. All individual 
weight maps are resampled and coadded in exactly the 
same way as the respective science images, yielding the 
noise properties for all pixels in the final coadded image 
(see Fig. for an example). 

4.6. The coaddition process 

Before the coaddition, all images are sky subtracted. In or- 
der to model the sky background we detect all objects in 
the field with SExtractor and replace them with the mean 
background as determined from the remaining pixels. This 
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Fig. 5. Coadded weight image of a small WFI@2.2 data 
set consisting of five exposures. One clearly identifies re- 
gions with less effective exposure time due to gaps be- 
tween CCDs and different pixel sensitivities. The size of 
the dither pattern also becomes obvious. Brighter regions 
correspond to pixels with higher weight. The variations 
from chip to chip are due to differences in the gain and 
the flatfield. 



image is then convolved with a very broad smoothing ker- 
nel (width between 200 and 500 pixels) and subtracted 
from the science image itself. 

For the coaddition the EIS drizzle in IRAF is used. It 
allows for a weighted mean coaddition, guaranteeing the 
best S/N in the stacked image. The resampling strategy 
ensures that the PSF is not artificially bloated in the stack, 
and that a varying pixel scale is correctly taken into ac- 
count also from a photometric point of view. Alternatively, 
one can use swarp (E. Bertin, Terapix), which contrary 
to drizzle, makes use of much more advanced resampling 
algorithms and does not lead to correlated noise in the 
stacked image. However, the differences between swarp 
and drizzle vanish if many images are stacked, which is 
the case in GaBoDS. 

Four factors determine the value of an output pixel 
in the coadded image. We have the input pixel value Jj 
from the science chip and an associated value Wi in the 
weight map. U represents the part of the input pixel that 
is mapped onto the corresponding output pixel. Besides, 
Jj is scaled with factors to the consistent photometric 
zeropoint and to a fixed exposure time (we chose 1 sec for 
this purpose): 

f . = 10 -0AZP i/ti) (3) 



where ti is the exposure time and ZPi the relative photo- 
metric zeropoint. All images are also weighted according 
to their sky noise. This weight scale is given by: 

w i = n 2 1 f 2 ■ ( 4 ) 

"sky.i^i 

Here we take into account the fact that the noise also 
scales with the flux scale /j. The values / out and W out in 
a stack of N images then read 

/out = ^= ff , W out = ^ WiWi . (5) 

EIS drizzle creates its output with the TAN projec- 
tion. Alternatively, th e COE projection can be used (see 
iGreisen fe CalabrettaL I2002L for further information on 
sky projections). In the stacked images North is up and 
East to the left. A reference coordinate can be specified for 
the coaddition. Thus, if multicolour information is avail- 
able for a particular pointing, the stacked images in the 
different bands are automatically registered with subpixel 
accuracy. 

5. Galaxy clusters in the background of NGC 300 

5.1. Characteristics of the data 

One of the ASTRO VIRTEL fields was a deep multi- 
colour observation of NGC 300 (see Table QJ, a face- 
on spiral galax y in the Sculptor grou p at a distance of 
about 2.1 Mpc l|Freedman et alll200l|) . Its angular size is 
25' x 18' , occupying about 40% of the WFI@2.2 field of 
view (34' x 34' ). The field did thus not meet all require- 
ments for GaBoDS, but the image seeing of a significant 
fraction of the i?-band data was around l'.'O. Deep V^-band 
observations (~ 37 ksec, 1'.'13 image seeing) were available 
too. This allowed us to search for distant galaxy clusters 
in the background of NGC 300 by means of weak lensing 
and the red cluster sequence method. 

Upon visual inspection in the Digitized Sky Survey 
(DSS) before the data retrieval request, we recognized 
two concentrations of fainter galaxies north-east and 
south-east of NGC 300, at the edge of the WFI@2.2 
fiel d. 21 s pectra were tak en for the first concentration 
bv lCappil jHeld fc Maranoj) (hereafter CHM98), confirm- 
ing a cluster at redshift 0.165. The second, less promi- 
nent concentrat ion is known as E DCC-499 at a red- 
shift of 0.117 (|Collins et all Il995[) . Hereafter we re- 
fer to these clusters as CL0056.03-37.55 (CL0053-37 in 
CHM98) and CL0056.02-37.90, respectively, using their 
epoch 2000.0 equatorial coordinates. The WFI@2.2 data 
for NGC 300 was taken for the ARAUCARIA project 
llPietrzvnski et al ., 2002a). an attempt to fine-tune the 
distance ladder by comparing different distance indicators 
such as cepheids, blue supergiants, tip of the red giant 
branch, and planetary nebulae for various nearby galax- 
ies. Since NGC 300 was monitored in 34 nights between 
July 1999 and January 2000, the creation of deep multi- 
colour images was feasible. 
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Fig. 6. PSF anisotropies in the NGC 300 i?-band image. Upper left: before correction, upper right: after correction. 
Lower left: Anisotropies as measured in the image, lower right: a two-dimensional polynomial fit to the PSF anisotropies 



5.2. Data reduction and catalogue creation 

The NGC 300 data was reduced essentially in the same 
way as described in Sect. 4. Due to the large extent of this 
galaxy a superfiat could only be calculated for the field 
outside NGC 300. Pixels lying inside the galaxy were only 
corrected for the gain differences, which were determined 
from the unaffected outer area. The same held for the sky 
subtraction of individual images before the coaddition pro- 
cess. The sky was modeled outside NGC 300 and assumed 
to be constant inside, so that no discrete jumps appeared 



between the inner and outer part. This image was then 
smoothed with a large kernel and subtracted. A substan- 
tial part of the images suffered from secondary scattering 
light, and from occasional vignetting caused by the filter 
holder. Most affected was the south-eastern corner of the 
field, where CL0056.02 is located. Our analysis therefore 
mainly concentrates on CL0056.03. 

The photometric zeropoint of the coadded U-band im- 
age was determined by matching stellar magnitudes to the 
seconda ry standard stars established bv lPietrzvnski et alJ 
in the Johnson-Cousins system. No calibration 
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was available for the i?-band, for which a zeropoint was 
determined based on the expected V — R colour for the 
red sequence of CL0056.03 (z = 0.165). We estimate it to 
be accurate within 0™1, which is sufficient for the analy- 
sis presented here since it does not rely on highly accurate 
absolute photometry. 

For the production of colour catalogues the coadded 
images were normalised so that their mean background 
noise erback = 1. They were then coadded without further 
rescaling, yielding a high signal-to-noise detection image. 
For this stack an adapted weight map was created com- 
bining the individual weight images accordingly. For the 
creation of a colour catalogue we used the detection im- 
age, the detection weight map and the unnormalised im- 
ages from filters V and R, as an input for SExtractor. Thus 
the flux for objects in different filters was measured within 
the same aperture, yielding relative colour information for 
galaxies with good internal accuracy. A further advantage 
of this approach was that the objects detected in the sin- 
gle catalogues for the V- and i?-bands were already in 
the same order and easily merged into one colour cata- 
logue. We calculated the galaxy colours from isophotal 
magnitudes (MAGJSO), and used MAG.AUTO for the 
magnitudes themselves. 

5.3. Using weak shear for cluster detection 

In the following standard weak lensing notations 
are used. For a technical review of this topic see 



iBartelmann fc Schneider! l)200l[) . The tidal gravitational 
field of a cluster-sized mass concentration induces a coher- 
ent distortion pattern in the images of distant background 
galaxies. By scanning the field for such characteristic dis- 
tortion patterns one can detect mass conc entrations di- 
rectly , irrespective of their luminosity (see |ErJjen el_alJ, 
200 0t lUmetsu fe Futamasel l200(j IWittman et all l200l[ 



Mirall es et all 120021 for examples). 



W e use the aperture mass statistics M ap l|Schneiderl 
1996) for the detection of galaxy clusters. M ap is a filtered 
integral of the projected mass distribution, k, inside an 
aperture. Its definition reads 



M a p = / d 2 tfK(tf,Z d ,Z s )£/(l?) . 
JO 



(6) 



Zd and z s are the lens and source redshifts, and t/(i?) is a 
compensated filter, i.e. J Q d-dflU (•$) = 0. The filter func- 
tion U({f) is chosen as 




(7) 



By defining the new filter function 



(8) 



M ap can be expressed in terms of the tangential shear j t , 
for which the observable ellipticities e t of the background 
galaxies are an unbiased estimator, as 



M, 



f d 2 tf 7t (tf)Q(tf) « - V e t (0i)Q(pi) • 
Jo n „■ 



(9) 



Hence, we can calculate the scalar M ap directly from ob 



ap 

servables. The noise for M ap is evaluated as 



J Map 



(10) 



where <t c is the ellipticity dispersion and n the num- 
ber density of background galaxies. We determined n = 
20 arcmhr 2 and a t = 0.34 for the WFI@2.2 data at hand. 
Based on this expression, the S/N ratio for mass peaks can 
be estimated. For doing so, a projected mass distribution 
Zd, Zs) at a redshift 2d is assumed, together with a 
redshift distribution p(z s ) for the lensed galaxies. The ex- 
pected signal S then reads 



S = 



dz s z d , z s )U(-d)p(z s ). 



(11) 



For t he mass d istribution a NFW profile ijNavarro et all 
Il996|) was used. lKruse fc Schneider! l)l999|) . showed how to 
convert this profile into the projected surface mass den- 
sity, k. For the source redshifts we took the normalised 
distribution 



2z \z 



exp 



1.5 



(12) 



proposed bv lBrainerd et alJ l|l996|) . We fixed the parame- 
ter zq based on photometric redshifts, computed from the 
UBVRI WFI02.2 data (EIS Deep Public Survey) for the 
Chandra Deep Field South. The exposure time in the R 
filter was 9000 seconds, lower than the ~ 15000 seconds 
for the NGC 300 field. Fig. □ shows the photometric red- 
shifts and the fit to the p(z s ) profile. A value of zq = 0.37 
is determined from the fit. We conservatively increased 
this value to zq = 0.4 for the S/N computation, given the 
significantly larger exposure time in R. 

Fig. [S] shows some predicted S/N ratios for massive 
haloes as a function of redshift. From this plot we see that 
clusters with masses of 3.2 x 10 14 A/ Q are detectable up to 
a redshift of z = 0.35 (0.42), and those with 1.0 x 10 15 M Q 
can be found up to z = 0.5 (0.65) for a filter size of 3'.2 
(4'.0, not shown). We are not sensitive to structures with 
~ 1.0 x 10 14 M Q at the 3'.2 scale, but can detect them at 
the 4'.0 scale up to z = 0.14. For comparison, S/N ratios 
are also shown for the case when background galaxies are 
at somewhat higher redshift (zq = 0.5). 

The S /N ratio for M ap peaks in a given data field can 
also be determined from the data itself, by randomizing 
ellipticities while keeping galaxy positions fixed. In the 
following we make use of this fact. For the weak lensing 
analysis the i?-band image was used. To avoid any bias- 
ing of the detection algorithm, NGC 300 was masked and 
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Fig. 7. Photometric redshift distribution for the Chandra 
Deep Field South, determined from UBVRI WFI@2.2 
photometry. The redshifts were es timated with hyperz 
(solid line) ijBolzonella et all l2000|) . Only objects that 
had a good redshift fit of P(x 2 ) > 0-9 went into the 
distribution shown. The dash-dotted line shows an inde- 
pendent measurement of the photometric redshifts from 
E. Hatziminaoglou (ESO) for the same data set (private 
communication). The peak at z — 2.8 is due to a de- 
generacy between low and high redshift galaxies and ap- 
peared since no prior for the luminosity of galaxies was 
used. These excess objects are randomly drawn from the 
low-z regime. Thus the main peak of this distribution 
was lowered, but not shifted in redshift. Equation (|12|l 
was fitted to both distributions, yielding consistent values 
z = 0.373 (0.377). 

replaced by the mean sky background, removing most of 
the flux present in the image. With SExtractor all ob- 
jects with at least 6 connected pixels > 2a above the 
sky background noise were detected. This catalogue con- 
tained 44146 objects. The shear estimates for these ob- 
jects wer e determined with the KSB algorithm as pro- 
posed bv lKaiser et alJ 1 19951) . An extensive description of 
our approach using the KS B method, includin g the PSF 
corrections, can be found in lErben et al.l l|200lj) . Figure 
shows the measured and the corrected PSF anisotropies 
for this particular field. After filtering, the background 
galaxy catalogue contained 12694 objects with securely 
determined shapes, a detection significance v > 12 (calcu- 
lated by KSB) and a bright magnitude limit of R > 23.0. 
This translates to the above mentioned number density of 
~ 20 galaxies per square arcmin. 

Fig. |U shows the M ap results for various filter widths. 
As can be seen, we recovered CL0056.03 in all but the 
largest filter scale at a 3cr level within 50-100" south of 



Fig. 8. The expected S/N ratios for various massive dark 
matter haloes, as detected with M ap in the NGC 300 
WFI@2.2 field. The effect of a background galaxy pop- 
ulation with higher redshift is shown, too. 



the cluster centre, whereas CL0056.02 was not detected. 
For the latter cluster useful shear information could only 
be obtained from about 30% of the area which would be 
available if the cluster did not lie next to the field corner 
in a region with bad image quality. This increased the 
noise in M ap by a factor of ~ 3. The argument of field 
truncation also holds for CL0056.03, but to a much lesser 
extent. There the useable field was limited by the edge of 
the image, 2'.5 east of the cluster centre, and NGC 300 
5'5 south-west. From our previous S/N considerations we 
conclude that a mass of ~ 1.4 x 10 14 M Q at redshift of 
z = 0.165 would produce a comparably significant Af ap 
detection. 

The M ap statistics furthermore picks up a number of 
other peaks. In the upper left panel of Fig. [U] the one 
for the smallest filter scale, we find a 3a peak that drops 
to 2.5cr in the upper right panel, and vanishes for larger 
filter scales (clump 'C'). Another peak below CL0056.03 
is detected at the 3a level for all but the largest filter scale 
(clump 'A'). Both peaks lie within 50" and 20-70" of two 
clumps of red galaxies, as is shown in the next section. 
There we determine their redshifts 2d, A = 0.47 and Zd,c = 
0.43. Our S/N calculations then yield masses of Ma(4'.0) rj 
(4 ± 3) x 10 14 M Q and M c (3'.2) w (6 ± 3) x 1O 14 M for 
these two clumps and for the filter scales in which they 
are detected most significantly. See also Fig. ^| for those 
two concentrations of red galaxies. 

Besides, a stable peak is found inside the mask of NGC 
300, with a detection limit of 3 — 3.5a for the three smaller 
filter scales. Such a detection is not surprising, since M ap 
is a highly non-local measure. It can pick up those parts 
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Fig. 9. The Af ap statistics for various filter scales: 3'.2 (upper left), 4'.0 (upper right), 5'.6 (lower left) and 8'.0 (lower 
right). Shown are M ap probability contours drawn from 5000 randomisations each. The contours depict the 2.0, 2.5, 
3.0, 3.5ct levels. The cluster CL0056.03 can be found left of the brightest star, and CL0056.02 sits in the lower left 
corner. See also Figs. El an d El f° r comparison. The large elliptical contour arises from the fact that M ap can not 
determine any value for the pixels inside the contour, since for those no galaxies lie inside the filter. Thus, the distance 
from this ellipse to the outer edge of the mask depicts the radius of the M ap filter function. 



of the shear field from a possible cluster hidden behind 
NGC 300 that extend beyond NGC 300 itself. Shear fields 
for massive clusters of g alaxies have been traced beyond 
10' of cluster centres fsee lClowe fc Schneiderll2002l for ex- 
amples), thus a sufficiently massive cluster behind NGC 
300 could easily be picked up. The detection in question 
is at the very outer edge of the galaxy disk, but the con- 



fusion limit of foreground stars in NGC 300 is already 
reached. Yet the optical thickness of the disk is still small, 
so that three larger and brighter isolated galaxies can be 
seen within a 2' wide window. Thus, if there was a mas- 
sive lower redshift cluster such as CL0056.03 at this po- 
sition, it could be seen through the disk. There is no in- 
dication for such an object at this position. More distant 
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Fig. 10. The (V — R, i?)-colour-magnitude diagram for all 
galaxies around NGC 300. The selected cluster sequence 
for CL0056.03 and CL0056.02 is marked by the large box. 
We do see clustering of redder objects in the upper window 
also. See the text for more details. 



Fig. 11. Shown are predictions for the WFI@2.2 V-R 
colours of elliptical galaxies as a function of redshift 
llBolzonella et all l200d: feruzual fc Charlotlll993|) . Based 
on the track that includes evolutionary effects we estimate 
the redshifts for the shear-selected red clumps. 



clusters, however, could no longer be identified as such, 
since their smaller and fainter images are drowned in the 
foreground confusion. No exces s X-Ray emission is found 
in the ROSAT All Sky Survey llVoges et allll999h for the 
M ap detections, apart from CL0056.03. 

5.4. Using the red sequence for cluster detection 

In a (V— R, R) colour- magnitude diagram fFig. lTUI) . plot- 
ted for all galaxies in the vicinity of NGC 300, we loosely 
selected a red cluster sequence with 0.67 < V — R < 
0.87,13.0 < R < 21.0. The sky distribution for objects 
inside this window is plotted in Fig. El (big red dots), 
whereas all other objects with R < 21 and outside the win- 
dow are shown as small dots. We then calculated the rms 
of the projected density of all red sequence members, and 
overlaid isodensity sigma-contours for their distribution. 
The smoothing length for the density was 3'.6. We recover 
CL0056.03 and CL0056.02 at the 8er and 5a level in over- 
density, respectively. Thus the cluster sequence we selected 
in the colour-magnitude diagram is actually a merger of 
two clusters of galaxies. 

The remaining galaxies that lied within the red se- 
quence window, but which did not belong to either of the 
two clusters, did not clump throughout the rest of the 
field on a level higher than la. Furthermore, all galaxies 
brighter than R — 20 and outside the red sequence win- 



dow showed no clumping at a level higher than 2a. We ad- 
ditionally checked the clustering properties of all objects 
brighter than R = 23, and did not find any highly signif- 
icant clumping apart from galaxies with 1.37 < V — R < 
1.56, 19.0 < R < 22.5. There we found four significant 
overdensities 'A' to 'D', with 'A' and 'C coincident with 
M ap detections (see above). Table |3] summarises the prop- 
erties and redshift estimates based on the track for ellipti- 
cal galaxies in Fig. We see that all clumps lie at simi- 



Table 3. All clumps except 'C lie close to the peak of the 
track in Fig. ^2 reducing the ambiguity in their redshift 
estimate. 



Red clump 


(V-R) 




A 


1.46 ±0.02 


0.47 ± 0.05 


B 


1.50 ±0.03 


0.47 ±0.05 


C 


1.41 ±0.02 


0.43 ±0.05, 0.51 ±0.05 


D 


1.50 ±0.03 


0.47 ±0.05 



lar redshifts, with some ambiguity for clump 'C. Actually, 
peaks 'A' and 'B' could form a connected system if spec- 
troscopically confirmed to be at the same redshift. More 
details are given in the caption of Fig. ED 
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Fig. 12. Galaxies inside the lower red cluster sequence of Fig.EHare shown as big dots. The small dots indicate galaxies 
brighter than R = 21 which do not fall inside the red sequence window. The overlaid contours are isodensity contours 
for the red sequence members, smoothed at a 3'.6 scale, and starting with the lc-overdensity contour in steps of la. 
CL0056.03 and CL0056.02 are detected at the 8a and 5a level respectively. CL0056.03 appears strongly elongated, 
with an intersecting 12' long filament extending north-south at its eastern side. Part of the filament could belong to 
Abell S0102 at (a,6) — (13.91,-37.41), a poor cluster at z = 0.05. Galaxies with measured spectra are highlighted 
with small squares around them. 



5.5. Butcher-Oemler effect in CL0056.03 

Fig. ^] depicts the two colour-magnitude diagrams for 
CL0056.03 and CL0056.02. Shown are galaxies with V - 
R = 0.80±0.5, 13.0 < R < 21.0 and within 5'.5 of the clus- 
ter centres. These galaxies were split into three distance 



bins (see Fig. 114(1 . in order to detect a possible evolution 
in their colour as a function of distance from the cluster 
core. The number of interlopers is small; we counted 41 
galaxies with the same properties in a 6'.5 x 33' strip on 
the opposite side of NGC 300, corresponding to a number 
density of about 0.19 galaxies per square arcmin. If no 
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further clumping takes place along the line-of-sight, 11% 
(19%) of interlopers can be expected within the area under 
consideration for CL0056.03 (CL0056.02). 

For CL0056.03, in the innermost bin 11 out of 19 galax- 
ies have measured redshifts, in the second bin it is 10 out 
of 60 galaxies, and in the outermost bin none of the 34 
galaxies have known redshifts. 10 galaxies appear to lie 
above the upper envelope (V — R > 0.87) of the clus- 
ter sequence, 7 of those are from the second distance bin, 
and the rest comes from the third bin. Their number is 
in agreement with the expected amount of interlopers. 4 
of those 10 galaxies are spatially concentrated within 1'.2. 
None of the 10 has spectroscopically determined redshifts, 
we assume that they are in the background of CL0056.03 
since they lie beyond the upper envelope of V — R < 0.87 
as defined by the red sequence. 

A closer inspection of the colour-magnitude diagram 
for CL0056.03 reveals several features: 

— Galaxies within 172 kpc of the cluster core (red dots) 
follow a tight correlation of 

V - R = (-0.013 ± 0.004)i? + 1.028 ± 0.077 (13) 

— Galaxies in the range [172-635] kpc of the cluster core 
(green triangles) show a steeper correlation with sub- 
stantially larger scatter: 

V - R = (-0.030 ± 0.014)i? + 1.376 ± 0.264 . (14) 

Less luminous galaxies with R > 19 are disproportion- 
ately bluer than those with R < 19, with the latter 
ones lying on the same sequence as the galaxies inside 
the innermost annulus. Fainter galaxies inside 172 kpc 
are bluer too, but to a much smaller extent than the 
galaxies in the second distance bin. 

— Galaxies within 635 kpc and R < 19 have V — R = 
0.795 ± 0.02, showing no evolution. 

— Galaxies in the range [635-927] kpc show no obvious 
concentration in colour-magnitude space, but the red 
sequence appears as an upper limit for the colour of 
those galaxies. Only two of them show significantly 
redder colours, whereas about a dozen appear bluer 
and about 20 lie on the red sequence as defined by 
galaxies in the inner two annuli. The upper limit in- 
dicates that a significant fraction of galaxies further 
away than 635 kpc could belong to the cluster popu- 
lation. However, without any measured redshifts this 
fraction can not be quantified at this stage. 

To summarise, for CL0056.03 a clear colour evolution for 
R > 19 and inside 635 kpc is found, increasing with the 
distance from the cluster core. Furthermore there exists a 
population of galaxies beyond 635 kpc which appears to 
share the same evolution limit as galaxies further inside, 
but they do not form a clear cluster sequence. The fraction 
of blue galaxies increases towards larger radii. 

For CL0056.02 we find similar, but less significant re- 
sults. The cluster sequence is broader and not so well- 
defined than the one for CL0056.03. There is no clear cut- 
off towards redder colours. Only galaxies inside the inner- 
most distance bin and with R < 19 form a tight sequence 



Table 4. The Butcher-Oemler Effect. In the case of 
CL0056.03 a galaxy is defined to be 'blue' if V-R < 0.755, 
i.e., if it lies more than 2a below the cluster sequence 

V - R = 0.795 ± 0.02, as defined by galaxies with R < 19. 
Galaxies with V — R > 0.87 are exluded from the statis- 
tics. For CL0056.02 the 2cr-threshold for 'blue' galaxies is 

V — R < 0.79. Galaxies with V - R > 0.85 are excluded 
from the statistics. Note that the number of interlopers is 
reduced in the outer two annuli since they are significantly 
truncated by the edge of the field of view. 



Distance from 


Galaxies 


Fraction of 


est. No. of 


cluster centre 


with(out) z 


blue galaxies 


interlopers 


CL0056.03 




0-200 kpc 


11 (8) 


0.26 


0.6 


200-740 kpc 


10 (50) 


0.35 


6.8 


740-1080 kpc 


(34) 


0.57 


6.5 


CL0056.02 




0-143 kpc 


1 (6) 


0.29 


0.6 


143-430 kpc 


1 (26) 


0.59 


4.5 


430-780 kpc 


(23) 


0.83 


6.2 



with V— R = 0.831±0.009. Galaxies in the second bin have 
a considerably larger scatter for brighter magnitudes, but 
a smaller one at the faint end as compared to CL0056.03. 
Objects within 385 kpc show a clear colour evolution with 
a slope of -0.029 ± 0.019 for R > 18.5. Beyond 385 kpc 
the galaxy population becomes significantly bluer, 83% 
of them lie below the red sequence. The fraction of blue 
galaxies is larger than for CL0056.03 in the outer two an- 
nuli, whereas in the cluster cores their fraction is compa- 
rable. We note that the data quality in the field around 
CL0056.02 is significantly lower than for CL0056.03, in- 
creasing the photometric errors. Table 01 summarises the 
Butcher-Oemler effect found for both clusters. 

5.6. Substructure in CL0056.03 

In CHM98 a velocity dispersion of l^l^kms -1 is 
stated for CL0056.03, based on the redshifts in Table 
From this dispersion the authors calculated a mass of 
~2x 1O 15 M0, assuming a virialised cluster and spherical 
symmetry. This value places CL0056.03 in the same league 
as supermassive clusters such as Abell 1689. Looking at 
Fig. El one notices the large elongation (~ 1 : 4) in the 
distribution of the red sequence members, arguing against 
spherical symmetry and relaxation, in which case the virial 
theorem yields inaccurate results. Given the large angu- 
lar separation of 105" (300 kpc) of the two cD galaxies, 
we searched for a correlation between the positions of 
galaxies and their redshifts, indicating an ongoing merg- 
ing process. The galaxies were split into two samples by 
drawing the perpendicular bisector of the connection line 
between the two cD galaxies (see left panel in Fig. I15|) . 
This way 10 galaxies with redshifts were assigned to the 
eastern cD galaxy (hereafter cDi), and 11 to the west- 
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Table 5. Heliocentric redshifts for CL0056.03, taken from 
CHM98. Here, v denotes the radial velocity in kms" 1 to- 
gether with its measurement errors a v , and the index m in- 
dicates whether a galaxy is closer to the eastern cD galaxy 
(ID = 6), or to the western one (ID = 16). 



ID 


RA (2000.0) 


DEC (2000.0) 


V 




z 


m 


1 


00:56:10.23 


-37:34:07.9 


51570 


282 


0.17202 


1 


2 


00:56:09.22 


-37:33:59.5 


50799 


51 


0.16945 


1 


3 


00:56:07.07 


-37:33:49.4 


51628 


138 


0.17221 


1 


4 


00:56:09.71 


-37:33:39.9 


50709 


120 


0.16915 


1 


5 


00:56:01.49 


-37:33:31.6 


49315 


111 


0.16450 


1 


6 


00:56:06.41 


-37:33:18.0 


49754 


154 


0.16596 


cD 1 


7 


00:55:57.50 


-37:32:51.5 


46839 


38 


0.15624 


2 


8 


00:55:58.90 


-37:33:00.7 


48490 


50 


0.16175 


2 


9 


00:56:02.76 


-37:33:03.9 


49588 


102 


0.16541 


1 


10 


00:56:04.76 


-37:33:00.3 


50394 


117 


0.16810 


1 


11 


00:56:04.75 


-37:32:51.0 


51594 


84 


0.17210 


1 


12 


00:55:55.59 


-37:32:20.2 


48514 


145 


0.16183 


2 


13 


00:55:56.74 


-37:32:35.4 


49555 


123 


0.16530 


2 


14 


00:55:57.74 


-37:32:23.4 


48256 


101 


0.16096 


2 


15 


00:55:59.32 


-37:32:35.2 


49041 


145 


0.16358 


2 


16 


00:55:59.09 


-37:32:18.2 


48888 


117 


0.16307 


cD 2 


17 


00:56:01.32 


-37:32:41.9 


50168 


187 


0.16734 


2 


18 


00:56:02.45 


-37:32:18.0 


47381 


138 


0.15805 


2 


10 


00:56:02.41 


-37:32:05.7 


49097 


113 


0.16377 


2 


20 


00:55:56.18 


-37:31:48.3 


49439 


14 


0.16491 


2 


21 


00:56:09.54 


-37:31:49.0 


49863 


107 


0.16633 


1 






1 75 


I 

<z> = 0.1685 ± 0.0029 

I * • 


<z> = 0.1624 ± 


0.0032 _ 





1 70 
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Fig. 16. The redshifts of the galaxies in CL0056.03 as a 
function of distance from the geometrical cluster centre. 
There is a highly significant (~ 4cr) correlation between 
galaxy position and redshift, indicating that CL0056.03 
actually consists of two possibly merging subclumps. The 
slope of the linear fit is different from zero on the 4.1a 
level. Note that cD x has a peculiar motion of 780 kms -1 
with respect to its own clump. 



ern one (CD2). Fig. 1161 shows the redshifts of galaxies 
as a function of distance from the geometric cluster cen- 
tre, defined as the centre of the connection line between 



Table 6. Properties for the eastern and western clumps 
around cDi respectively cD 2 . 



Clump 


<*> 


(v) [kms -1 ] 


(?v [km s 1 ] 


cDi 


0.1685 ± 0.0029 


50521 


879 


cD 2 


0.1624 ± 0.0032 


48697 


960 



cDi and CD2. We find different mean redshifts for the 
galaxies around cDi and CD2, translating into a veloc- 
ity difference of ~ 1824 km s _1 . The velocity dispersions 
of the two clumps are thus significantly smaller than the 
one given in CHM98 for the whole system, as is the to- 
tal dynamical mass. Details can be found in Table [BJ In 
order to obtain a more rigorous estimate for the signif- 
icance of this discrepancy, we used Monte Carlo simu- 
lations. Galaxy redshifts were randomly distributed 10 6 
times while keeping galaxy positions fixed, destroying any 
possible correlation between redshift and position. In only 
95 out of the 10 6 cases do we find a higher mean redshift 
for the clump around cDi and a lower one for the other 
than in the observed data, meaning a 3.9<t significance. In 
addition, we checked whether the two distributions with 
&v,i — 879 km s -1 and a Vi % — 960 km s^ 1 can be drawn 
from the same parent distribution with a v = 1296 km s^ 1 . 
We created 10 5 random realisations of this parent gaus- 
sian distribution, each containing 21 velocities that were 
ordered. Every realisation was split into two sub-samples, 
having the same number of overlapping objects in velocity 
space as the real data, shown in Fig.^3 These sub-samples 
contained 10 and 11 mock galaxies, and it was checked 
whether the difference in their mean velocities was larger 
than the observed 1824 km s" 1 . We found that in 96.9% 
of all cases the observed velocity histogram could not be 
drawn from a single gaussian distribution, thus the clumps 
probably did not yet mix. 

We note that the velocity difference between cD\ and 
cD 2 themselves, however, is v — 866 km s -1 and thus 
much lower than the difference between the sub-clumps. 
We conclude that the 21 measured redshifts for CL0056.03 
are probably not representative for the ~ 113 cluster mem- 
ber candidates, as identified by their colours. In the fol- 
lowing a simple dynamical model is used in order to check 
whether the system is gravitationally bound or unbound. 
A linear orbit of the clumps is assumed, i.e. there is no 
shear or rotation component. The system is seen under 
some inclination angle 95, with ip = if the clumps were 
aligned along the line of sight. Based on Newtonian dy- 
namics the condition of a bound system can be written as 
v 2 r < 2GM, or, taking into account the inclination angle 
ip, as 

*5 p < sin 2 (j9 cos <p , (15) 
2GM y y ' v ' 

where i> r = 1824 km s^ 1 is the radial velocity difference 
between the clumps and r p = 300 kpc the projected dis- 
tance on the sky between cDi and CD2. M is the total 
mass of the system. The right hand side of (fT5|l can not 
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get larger than 0.385 for ip = 54°. For M < 3.0 x 10 14 M 
the expression ((ToT) does not allow bound solutions for the 
given v r and r p . For M — 5.0 x 10 14 M the probability for 
a bound system is 49% (the fraction of allowed inclination 
angles over all inclination angles) , and it increases to 70% 
for M = 1.0 x 10 15 M©. Thus, the lensing estimated mass 
of M = 1.4 x 10 14 M© for CL0056.03 is not sufficient for a 
bound solution. If the velocity difference v r = 866 km s -1 
between the two most luminous galaxies cDi and CD2 was 
representative for CL0056.03, instead, then in 57% of the 
cases the system would be gravitationally bound. 

Our interpretation for the dynamical state of 
CL0056.03 is thus that the eastern clump is closer to us 
and the western one further away, and that the system is 
in a pre-collision phase with the two clumps approaching 
each other, since probably no mixing has yet taken place. 
In addition, there is some indication from the imaging data 
itself, arguing for this scenario: both the eastern and the 
western clump are embedded in large, but well separated 
haloes, each measuring about 340 kpc. The halo around 
the western clump, in particular, shows tidal features. We 
take this as an evidence for ongoing tidal stripping and 
merging processes inside the clumps, but probably no en- 
counter has yet taken place between them. In the spectra 
taken by CHM98 no emission lines are found for the clus- 
ter members, thus there is no sign for nuclear activity 
or major star formation. The other scenario is that there 
was already an encounter and the two clumps are reced- 
ing, which would place the western clump closer to the 
observer. Based on the small number of spectra the latter 
state can not be ruled out. 

Besides the two sub-clumps there is some evidence for 
a more extended structure in CL0056.03. In Fig. ^]a fil- 
ament is seen, extending to about 7' north and 4' south 
of the cluster centre, but not lined up with the cluster 
centre itself. The angle between this linear filament and 
the orientation of CL0056.03 is ~ 40 degrees. North of 
CL0056.03 at (a = 00 : 55 : 45.6, 5 = -37 : 24 : 46), how- 
ever, the small cluster Abell S0102 is found at a redshift 
of z = 0.056. A contamination of the possible filament 
by member galaxies of Abell S0102 can not be ruled out. 
Furthermore, 5'.2 from the cluster centre and in the direc- 
tion of the cluster orientation, an isolated elliptical galaxy 
with R = 14.4 and V - R = 0.745 is found, 0.05 mag 
bluer than the cluster sequence. Its distance to the centre 
of Abell S0102 is 6'.2. Based on the present data it can not 
be concluded whether this elliptical is isolated in space or 
a physical member at the fringes of one of the two clusters. 

6. Conclusions 

We introduced the Garching-Bonn Deep Survey, a 12 
square degree survey for weak lensing investigations. It 
is primarily a 'virtual' survey, since about 75% of the 
data was taken from the ESO archive. Within a dedicated 
ASTRO VIRTEL program the functionality of the querator 
search engine was substantially expanded. It now allows 
an effective filtering of the data with respect to charac- 



teristic properties. One of these archival fields, centreed 
on NGC 300, showed two concentrations of background 
galaxies upon visual inspection of the field in the DSS. 
The cluster nature of these two concentrations was spec- 
troscopically confirmed beforehand by other groups. 

Making use of the high quality of the i?-band expo- 
sures, we showed that one of the two clusters (CL0056.03) 
can be detected by its weak gravitational shear signal. 
The second cluster (CL0056.02) was not found in the weak 
lensing data, probably due to field truncation and insuf- 
ficient image quality. Besides, two other significant and 
coherent shear patterns were detected. Both of them coin- 
cide with significant overdensities of red galaxies, for which 
we gave redshift and mass estimates, based on their appar- 
ent colour and lensing strength. A third shear detection 
was seen behind NGC 300 itself, but there is no further 
evidence from background light for a cluster of galaxies at 
this position. Thus, apart from this hidden detection, all 
Map peaks equal to or higher than 3er in filter scales > 
3'.2 coincide with overdensities of red galaxies. These M ap 
peaks, however, are at the limit of what can be concluded 
from this data set. The tidal gravitational fields of more 
massive clusters of galaxies imprint shear fields in the im- 
ages of background galax ies that are easily detecte d on 
the 10cr level and beyond l|Clowe fc SchneiderL I2OO2I . 

Based on the high virial mass of CL0056.03, derived 
from spectroscopic data by CHM98, and the compara- 
tively shallow weak lensing signal, a closer analysis of this 
cluster was performed. A highly significant correlation be- 
tween redshift and the position of galaxies in the clus- 
ter is found, indicating the existence of two sub-clumps, 
separated by 1824 km s -1 in velocity space. Based on the 
present data and the derived lensing mass it can not be 
concluded whether the system is gravitationally bound. 
By looking at the projected sky distribution of galaxies 
with similar properties in a (V — R, R) colour-magnitude 
diagram, we find an extended filament of galaxies inclined 
to the cluster's major axis. Without further spectroscopic 
data of galaxies in this filament its nature can not be fur- 
ther clearified. We take the presented properties as evi- 
dence that CL0056.03 is a younger cluster in formation 
and has not yet reached its equilibrium state. 

CL0056.03 shows a clear and tight red sequence in 
the (V — R, R) colour-magnitude space, with an increas- 
ing fraction of bluer galaxies towards larger cluster radii. 
CL0056.02 also shows a red sequence, the results presented 
for this cluster, however, are more uncertain due to field 
truncation and inferior data quality. 
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Redder objects in the NGC 300 field 




14.0 13.8 13.6 13.4 

RA [degree] 



Fig. 13. Shown is the distribution of galaxies inside the upper red cluster sequence of Fig. EH (big dots). The small 
dots indicate galaxies brighter than R — 22.5 which do not fall inside the red sequence window. Clump 'A' is a very 
tight concentration of nine galaxies with R < 22.3 within 25" . There are five more galaxies spread along a 1' filament 
to the east and to the west. Clump 'B' is a much looser, 4' long filament consisting of 13 galaxies without a central 
concentration. In the centre of the filament one finds a R = 19.4 bright elliptical, about 1 magnitude brighter than 
the second brightest member of this concentration. Clump 'C consists of 6 galaxies within 80" , the brightest one 
with a R-magnitude of 20.9. Finally, clump 'D' is a very loose grouping of fainter red galaxies that looks like a chance 
alignment rather than a cluster. Note that the resolution of this plot is not high enough to reveal all candidate galaxies. 
The smoothing length for the density contours was 1'.5 arcmin. Clumps 'A' and 'C lie within ~ 50" of weak lensing 
peaks. 



20 



M. Schirmer et al.: GaBoDS: Anatomy of galaxy clusters in the background of NGC 300 

CL0056. 03-37. 55 CL0056. 02-37. 90 



.2 - 



1 .0 - 



X A 



0.8 - EP H H 



0.6 - 

- • < r < 1 72 kpc 

0.4 - a 1 72 < r < 635 kpc 

- x 635 < r < 927 kpc 

(100 kpc = 35.4") 
0.2 \ 




14 15 1i 



17 18 
R 



20 21 



0.8 



0.6 



0.4 



0.2 




. < r < 128 kpc 

a 128 < r < 385 kpc 
x 385 < r < 698 kpc 

(100 kpc = 47.5") 



1 4 



15 



17 18 
R 



I 9 



20 21 



Fig. 14. Colour-magnitude diagrams for galaxies inside the cluster sequence window. Galaxies with spectroscopically 
determined redshifts are marked with squares, the cD galaxies are indicated with a double square (they also have 
measured redshifts). The galaxy with R = 19 and V — R = 0.73, marked with a diamond symbol, has z = 0.27, higher 
than the cluster redshift. Furthermore, galaxies were split into three distance bins as seen from the geometrical cluster 
centre. The outer radii (635 kpc respectively 385 kpc) of the second annuli were chosen in a way that galaxies in 
the third, outer annulus do not show an apparent concentration with respect to the red cluster sequence any more. 
Isodcnsity contours in the colour-magnitude space were calculated from all galaxies. 
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Fig. 15. Image of CL0056.03. In the left panel the 21 spectroscopically confirmed cluster members are shown, together 
with the perpendicular bisector of the connection line between the two cD galaxies, which we chose to split the cluster 
galaxies into two samples. An [Oil] emission line galaxy at higher redshift is indicated. The right panel shows a steeper 
scaled version of the image at left, showing that cDi and cD 2 are both embedded in large, but well separated haloes. 



